The Ou (3 /3-decay nuclear matrix elements with self-consistent short-range correlations 
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A self-consistent calculation of nuclear matrix elements of the neutrinoless double beta decays 



(OvP/3) of 7b Ge, S2 Se, Jh Zr, 



'Mo. 



116 Cd, 128 Te, 130 Te and 



Xe is presented in the framework 



of the renormalized quasiparticle random phase approximation (RQRPA) and the standard QRPA. 
The pairing and residual interactions as well as the two-nucleon short-range correlations are for the 
first time derived from the same modern realistic nucleon-nucleon potentials, namely from charge- 
dependent Bonn potential (CD-Bonn) and the Argonne V18 potential. In a comparison with the 
traditional approach of using the Miller-Spencer Jastrow correlations matrix elements for the Of/3/3 
decay are obtained, which are larger in magnitude. We analyze the differences among various 
two-nucleon correlations including those of the unitary correlation operator method (UCOM) and 
quantify the uncertainties in the calculated O^/3/3-decay matrix elements. 

PACS numbers: 23.10.-s; 21.60.-n; 23.40.Bw; 23. 40. He 

Keywords: Neutrino mass; Neutrinoless double beta decay; Nuclear matrix element; Quasiparticle random 
phase approximation 



o 

=5 
S3 



m 
O 

(N 
O 

o 



X 



I. INTRODUCTION 

The present data on neutrino oscillations shows that 
pattern of neutrino masses and mixing (Pontecorvo- 
Maki-Nakagava-Sakata mixing matrix) is not an analogue 
of that for quarks (Cabibbo-Kobayashi-Maskawa quark 
mixing matrix) [l], Q • The generation of neutrino masses 
can be explored, if the absolute scale of neutrino masses 
will be fixed and the issue of leptonic CP violation will be 
understood [|[ . This might happen, if the lepton number 
violating neutrinoless double beta decay (Of/3/3-decay) 
will be undoubtable observed in running [J, |5( or planned 
0, H, 0, HI Oi//3/3-decay experiments. 

The 0^/3/?-decay is a very sensitive probe for the Ma- 
jorana neutrino mass @, H, E3, EH The O^/3/3-decay 
can occur through different processes but all of them re- 
quire that the neutrino has nonzero mass and is a Ma- 
jorana particle [l3[. The most proximate or discussed 
theoretical model is to mediate the O^/3/3-decay by the ex- 
change of a light Majorana neutrinos. A measurement of 
the decay rate, when combined with neutrino oscillation 
data and a reliable calculation of nuclear matrix elements 
(NMEs), would yield insight into all three neutrino mass 
eigenstates, the type of neutrino mass spectrum (normal 
hierarchy or inverted hierarchy) and possibly Majorana 
CP-violating phases. 

An important subject in neutrino physics is a reliable 
calculation of the 0^-decay NME M 0v [3. Unfortu- 
nately, there are no observables that could be directly 
related to magnitudes of NMEs. The most popular nu- 
clear structure methods, which has been applied for this 
task, are proton-neutron Quasiparticle Random Phase 
Approximation (QRPA) with its variants [15| and the 
Large-Scale Shell Model (LSSM) [IE E3, El- Recently, 
there has been significant progress towards the reduc- 



tion of uncertainty in the calculated NMEs [HI, Ell . A 
detailed anatomy of the O^/3/3-decay NMEs pointed out 
a qualitative agreement between results of the QRPA- 
like and LSSM approaches It particular, it was 

shown that only internucleon distances r.y 2 — 3 fm 
contribute to M 0l/ , what explains a small spread of re- 
sults for different nuclei. Further, it has been shown that 
correlated NME uncertainties play an important role in 
the comparison of 0V/3/3-decay rates for different nuclei 

The improvement of the calculation of the nuclear ma- 
trix elements is a very important and challenging prob- 
lem. The problem of the two-nucleon Short-Range Cor- 
relations (SRC) have recently inspired new O^/3/3-decay 
studies Ho, HH, [H[. In the majority of previous calcu- 
lations SRC have been treated in a conventional way 
via Jastrow-type correlation function in the parametriza- 
tion of Miller and Spencer [24|. Recently, it was found 
that the consideration of the unitary correlation opera- 
tor method (UCOM) leads to increase of the O^/3/3-decay 
NME by about 20-30% [H E2, H- It was concluded 
that we do not know the best way to treat the SRC, a 
fact that contributes to uncertainties. 

In the present article we improve on the Miller-Spencer 
Jastrow and the UCOM and perform a selfconsistent cal- 
culation of the Oi^/3/3-decay NMEs by considering pair- 
ing, ground state and short-range correlations deduced 
from the same realistic nucleon-nucleon (NN) interac- 
tion. In particular, the two-nucleon short-range corre- 
lations will be determined within the coupled cluster or 
exponential S approach by using CD-Bonn and Argonne 
V18 NN-forces [26|, H3| and compared with Jastrow and 
UCOM SRC. Then, they will be used in the QRPA and 
renormalized QRPA (RQRPA) calculations of the 0i/(3/3- 
decay NMEs of experimental interest. 
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The paper is organized as follows: In Sec. II the for- 
malism of the 0z//3/3-decay associated with exchange of 
light Majorana neutrinos is briefly reviewed. Sec. II 
is devoted to the analysis of different treatments of the 
two-nucleon short-range correlations in the context of the 
correlated Oz^/3/3-decay operator. In Sec. IV we present 
numerical results for nuclei of experimental interest. Sec- 
tion V summarizes our findings. 



II. FORMALISM 

In this section we present basic expressions associated 
with the calculation of the Oz^/3/3-decay NME, what allow 
us to discuss the effects of Finite Nucleon Size (FNS) and 
two-nucleon SRC. 

By assuming the dominance of the light neutrino mix- 
ing mechanism the inverse value of the O^/3/3-decay half- 
life for a given isotope (A, Z) is given by 



f 



= G^(E o ,Z)\M ,Ov \ 2 \(m P0 )\ 



(1) 



Here, G 0v (Eq, Z) and M' u are, respectively, the known 
phase-space factor (Eq is the energy release) and the 
nuclear matrix element, which depends on the nuclear 
structure of the particular isotopes (A,Z), (A 7 Z + 1) 
and (A, Z + 2) under study. Under the assumption of the 
mixing of three light massive Majorana neutrinos the ef- 
fective Majorana neutrino mass (mpp) takes the form 



N 

E 



|C/ ei |V Ql m 4 , (all m, > 0) 



(2) 



where U e i is the first row of the neutrino mixing matrix 
and the and are unknown Majorana phases. It is ap- 
parent that any uncertainty in M' v makes the value of 
{rn-pp} equally uncertain. 

Our phase space factors G 0u (Eq, Z), which include 
fourth power of axial-coupling constant qa = 1.25, are 
tabulated in Ref. [281 ] ■ They agree quite closely with 
those given earlier in Ref. [29]. The G 0v (E ,Z) con- 
tain the inverse square of the nuclear radius (R n uci)~ 2 , 
compensated by the factor R nuc i in M' 0v . Different au- 
thors use different conventions for R n ud {Rnuci = roA 1 / 3 
ro = 1.2 fm or ro = 1.1 fm), a fact that is important 
to keep in mind when comparing the matrix elements 
without also looking at G°"(E , Z). 

The nuclear matrix element M' v is defined as 
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where M 0u consists of Fermi, Gamow- Teller and tensor 
parts as 



M F 

"A 



M G7 



M 7 



(4) 



This definition of M' v [15} allows to display the effects 
of uncertainties in 17^4 and to use the same phase factor 
G°"(Eq,Z) when calculating the Qi/(3(3-decay rate. 

In the QRPA (and RQRPA) M 0v is written as sums 
over the virtual intermediate states, labeled by their an- 
gular momentum and parity J n and indices ki and kf 
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The reduced matrix elements of the one-body operators 
c+c n (c„ denotes the time- reversed state) in the Eq. ([5]) 
depend on the BCS coefficients Ui,Vj and on the QRPA 
vectors X, Y [H[. The difference between QRPA and 
RQRPA resides in the way these reduced matrix elements 
are calculated. 

The two-body operators Ok, K — Fermi (F), Gamow- 
Teller (GT), and Tensor (T) in §5§ contain neutrino po- 
tentials and spin and isospin operators, and RPA energies 



F (r 12 ,E%) 
GT (r 12, E%) 
T (r 12 ,E%) 



r+(l)r+(2)H F (r 12 ,E%), 
r+(l)r+(2) H GT (r 12 ,E%) o 12 , 
t+(1)t+(2) H T (r l2 ,E%) S 12 - (6) 



with 
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S12 = 3((?i • r 12 )(a2 ■ ri 2 ) - oi 2 . 



(7) 



Here, ri and r 2 arc the coordinates of the nucleons un- 
dergoing beta decay. 

The neutrino potentials are integrals over the ex- 
changed momentum g, 
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The functions f F , G r(qr '12) = jo(qri 2 ) and fr(qr 12) = 
~32(q r i 2 ) ai 'e spherical Bessel functions. 

The potentials ([5]) depend explicitly, though rather 
weakly, on the energies of the virtual intermediate states, 
E k ]7r . The functions liK(q 2 ) that enter the Hk's through 
the integrals over q in Eq. ([8]) are 
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Here, we used the Partially Conserved Axial Current 
(PCAC) hypothesis. 

The FNS is taken into account via momentum depen- 
dence of the nucleon form-factors. For the vector, weak- 
magnetism and axial-vector form factors we adopt the 
usual dipole approximation as follows: 



9v(q 2 ) 

9M(q 2 ) 

9A(q 2 ) 



9v 



(l + q 2 /M 2 ) 2 ' 

(m p - Hn)gv(q 2 ), 

9A 
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(10) 



where gy = 1, = 1-00 (quenched) and 1.25 (un- 
quenched), (jj, p — //„) — 3.70. The parameters My = 
850 MeV and M A = 1 086 MeV come from electron 
scattering and neutrino charged-current scattering exper- 
iments. 

The 0^/3/?-decay matrix elements were usually calcu- 
lated in some approximations, which are only partially 
justified (see also discussion in (l5j): 
i) The effect of higher order terms of nucleon currents 
was not taken into account. In this case we have 



h F (q 2 ) = g v (q 2 ), 



(9 2 ) 



9 2 A (1 2 ) 
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(11) 

We note that if in addition nucleons are considered to be 
point-like hp and Hgt are equal to unity, 
ii) The closure approximation for intermediate nuclear 
states was considered by replacing energies of interme- 
diate states [Ej„ — (Ei + Ej)/2] by an average value 
I«10 MeV. 

Within these approximations the neutrino potential in 
Eq. ([5]) can be written as [30l | 



Hbare(r 12, E) 



- [sin (Er 12 )Ci(Er 12 ) 
cos (Er 12 )Si(Er 12 )] 



''12 



(12) 

Here, Ci(x) and Si{x) are the cosine and sine inte- 
grals, respectively. The value of E has practically no 
impact on the behavior of neutrino exchange potential 
at short internucleon distances. In the limit E = and 
zero neutrino mass the neutrino potential is coulombic: 

Hbare(ri 2 , E = 0) = R n ucl/fl 2 . 

It is worth to mention some general properties of the 
Fermi Mp and the Gamow- Teller Mqt matrix elements, 



in particular some multipole contributions of states of the 
intermediate odd-odd nucleus are equal to zero. We have 



M F (J H 
M F (J- 
M GT (0 A 



for odd J. 
for even J, 
0. 



(13) 



III. 



SHORT RANGE CORRELATIONS FOR 

THE O^/3/3-DECAY 



An important component of the Mr in © is an unan- 
tisymmetrized two-body matrix element, 



(p(l),p , (2); l /|| Ok II n(l), n'(2); J), 



(14) 



constructed from two one-body matrix elements by cou- 
pling pairs of protons and neutrons to angular momen- 
tum J . We note, in the closure approximation, i.e., if 
energies of intermediate states {Ej^ — Ei) are replaced 
by an average value E, and the sum over intermediate 
states is taken by closure, J2 n \ ^k ><: = 1> we en< ^ 
up with antisymmetrized two-body matrix elements. As 
the virtual neutrino has an average momentum of ~ 100 
MeV [20( , considerably larger than the differences in nu- 
clear excitation, the closure approximation limit is found 
to be meaningful showing on the importance of correla- 
tions of the two /^-decaying nucleons. 



The Jastrow and UCOM short-range 
correlations 



The QRPA (RQRPA) as well as the LSSM approaches 
do not allow to introduce short-range correlations into 
the two-nucleon relative wave function. The traditional 
way is to introduce an explicit Jastrow-type correlation 
function f(r\ 2 ) into the involved two-body transition ma- 
trix elements 



J 



firi 2 )0 K (r 12 )f(n 2 ) || tfy) 
|| K (r 12 ) || q/j). 



Here. 



I*7> 
l*7> 



/(n 2 ) 1*7), 

|n(l),n'(2);J) 



(15) 



(16) 



are the relative wave function with and without the short- 
range correlations, respectively. In the parametrization 
of Miller and Spencer [24| we have 

f(r 12 ) = l- e - ar2 (l-br 2 ), 0= l-lfm- 2 , b = 0.68 /m~ 2 . 

These two parameters (a and b) are correlated and chosen 
in the way that the norm of the relative wave function 
\^>j) is conserved. 

Usually, the nuclear matrix element M 0l/ is calculated 
in relative and center-of-mass coordinates as the Jastrow 
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correlation function depends only on r\2 . This is achieved 
with help of the well-known Talmi-Moshinski transforma- 
tion [3l[ for the harmonic oscillator basis. Within this 
procedure the chosen construction of the relative wave 
function, namely a product of /(V12) with harmonic os- 
cillator wave function in (|16p . is well justified. Any more 
complex structure of correlation function, e.g., a con- 
sideration of different correlation functions for different 
channels, would result in violation of requirements in Q13p 
as the Talmi-Moshinski transformation is considered. 

Recently, it was proposed [23] to adopt instead of the 
Jastrow method the UCOM approach for description of 
the two-body correlated wave function 25:]. This ap- 
proach describes not only short-range, but also central 
and tensor correlations explicitly by means of a unitary 
transformation. The state-independent short-range cor- 
relations are treated explicitly while long-range correla- 
tions should be described in a model space. Applied to 
a realistic NN interaction, the method produces a cor- 
related interaction, which can be used as a universal ef- 
fective interaction, for calculations within simple Hilbcrt 
spaces. The UCOM method produces good results for 
the bindin g en ergies of nuclei already at the Hartree- 
Fock level [32|]. There are also some first applications 
for description of collective multipole excitations [33T |. 

Within the UCOM the short-range and long-range 
correlations are imprinted into uncorrelated many body 
states by a unitary transformation. In the case of the 
O^/3/3-decay calculation the correlated two-nucleon wave 
function was taken as 

\*j) = Cr|*7>- (18) 

Here, C r is the unitary correlation operator describing 
the short-range correlations. The explicit form of C r is 
given in [25| with a separate parametrization for different 
LS-channels. In application to the O^/3/3-decay this fact 
leads to a slight violation of conditions (fl~3)) when Talmi- 
Moshinski transformations are considered. The UCOM- 
corrected NMEs of the Of /3/3-decay are significantly less 
suppressed when compared with those calculated with 
Jastrow SRC jMUl]. 



B. Self-consistent two-body short-range 
correlations 

The two-nucleon wave function with short range cor- 
relations can be calculated from the same realistic NN- 
interaction, which is used in the derivation of the Brueck- 
ner G-matrix elements of the nuclear Hamiltonian. A 
method of choice is, e.g., the Brueckner-Bethe hole- line 
expansion, the coupled cluster method (CCM) or expo- 
nential S approach and the approach of self-consistent 
evaluation of Green functions [261 ] . 

There are various modern NN potentials, which yield a 
very accurate fit to the nucleon-nucleon scattering phase 
shifts. Two of them are the so-called charge-dependent 



Bonn potential (CD-Bonn) [34[ and the Argonne V18 po- 
tential ( Argonne) [35[ . They differ in description of both 
long-range and short-range parts of NN-interaction. The 
CD-Bonn is derived in the framework of the relativistic 
meson field theory. The Argonne potential is a purely 
local potential in the sense that it uses the local form of 
the one-pion exchange potential for the long-range part 
and parametrizes the contributions of medium and short- 
range distances in terms of local functions multiplied by 
a set of spin-isospin operators. 

We have chosen the CCM [27| to evaluate the effect 
of short-range correlation on the Of/3/3-decay NMEs, be- 
cause it provides directly correlated two-body wave func- 
tions. The basic features of the CCM are described in the 
review article by Kiimmel (3?j] • The developments of this 
many-body approach with applications can be found in 

The CCM starts by assuming an appropriate Slater 
determinant \<fr) as a first approximation for the exact 
eigenstate of the A-particle system. The many-body 
wave function of the coupled cluster or exp(S) method 
can be written as 

\$) corr . = exp(^2S^J |$>. (19) 
The n-particle n-hole excitation operator S n is given by 



^2 (pi ■ ■■Pn\S n \i'i ■ ■ ■ v n )o) pi . . . a} pn a Vn ■ ■ -a vi . (20) 

The sum in (|20p is restricted to states pi which are un- 
occupied in the model state |$), while states z/j refer to 
states which are occupied in |$). 

A Slater determinant of harmonic oscillator wave func- 
tions is considered for |$). For a given nuclear system 
of interest an appropriate value of the oscillator length 
b is chosen. In the so-called S% approximation of the 
CCM the amplitudes defining S n with n > 3 in (jTSJ) 
are ignored. This means that effects beyond Hartree- 
Fock and two-body correlations (i.e. genuine three- and 
more-particle correlations) are ignored. This leads to a 
coupled set of equations for the evaluation of the correla- 
tion operators Si and §2 (27[. Therefore this S2 approxi- 
mation corresponds essentially to the Brueckner-Hartree- 
Fock (BHF) approximation of the hole-line expansion or 
Brueckner theory. In fact, the hole-hole scattering terms, 
which arc included in the S2 but ignored in BHF turn out 
to yield small effects only. Therefore it is consistent to 
combine the correlation effects from CCM with the ma- 
trix elements of the G-matrix, the effective interaction 
determined in the BHF approximation. 

The use of the oscillator ansatz in Slater determinant 
|$ > in cq. (Tr9T) leads to an evaluation of the correlated 




FIG. 1: Two nucleon wave functions as a function of the rel- 
ative distance for the 1 S'o partial wave and radial quantum 
numbers n = 0, 1, 2, 3 and 4. The results in panel a) are for 
uncorrelated two-nucleon wave functions. The results in pan- 
els b), c) and d), respectively, are with coupled cluster method 
with CD-Bonn potential, coupled cluster method with Ar- 
gonne potential and Miller-Spencer Jastrow short-range cor- 
relations. The harmonic oscillator parameter b is 2.18 fm. 



two-nucleon wave functions in terms of product wave 
functions for the relative and center-of-mass coordinates. 
The two-body states take the form 



[n(lS)j] NLJi 



(21) 



Here, N and L denote the harmonic-oscillator quantum 
numbers for the center of mass wave function and I refers 
to the orbital angular momentum for the relative motion, 
which is coupled with a total spin of the pair S to angular 
momentum J . The basis states for the radial part of this 
relative motion are labeled by a quantum number n. 

As an example we present in Fig. [T] relative wave 
functions for correlated and uncorrelated two-body wave 
functions in the case of 1 Sq partial waves and different 
values of the radial quantum numbers n (n = 0, 1, 2, 3 
and 4) . In panel a) uncorrelated harmonic oscillator wave 
functions are plotted. Panels b) and c) show the relative 
wave functions obtained with help of CCM employing 
the CD-Bonn and the Argonne potentials. For a com- 
parison relative wave functions with Miller-Spencer Jas- 
trow SRC are displayed in panel d). While the Jastrow 
ansatz suppresses the relative wave function in the limit 
fi2 completely, we find that this suppression effect 
is much weaker in the CCM calculation. This is true 
even if the Argonne potential is used, which is known 
to produce stronger short-range components than the 
softer CD-Bonn potential. Also note that the correlated 
wave functions derived from realistic interactions exhibit 
a short-range behavior which depends on the radial quan- 
tum number n, whereas the Jastrow approach yields al- 
most identical relative wave functions for small values of 
r\2- 
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FIG. 2: The ratio of neutrino potentials with and with- 
out two-nucleon short-range correlations (SRC). Results are 
shown for the CCM CD-Bonn and Argonne and Miller- 
Spencer SRC with and without consideration of the effect 
of finite size of a nucleon. It is assumed E — 8 MeV. 



components of the NN interaction at short distances 
are weaker for the CD-Bonn potential in comparison with 
Argonne interaction. But, in the case of the Jastrow 
SRC the reduction of the relative wave function for small 
values of r\i is even much stronger. 

The advantage of the CCM 27[ is a factorization of 
the correlated two-body wave function on a product of a 
correlation function and a harmonic oscillator wave func- 
tion. This allows us to discuss the effect of the SRC 
in terms of the correlated operator, which is a product 
of the transition operator Ok(t\i) and two correlation 
functions f(r\2) (see Eqs. (fT5|) V For our purposes we 
consider CCM CD-Bonn fs and CCM Argonne Ja corre- 
lation functions deduced from the ^c^n. = 0) correlated 
two-body wave function. The use of this single correla- 
tion function for all partial waves and quantum numbers 
n is numerically well justified and is dictated by the use 
of the Talmi-Moshinski transformation in the evaluation 
of the O^/3/3-decay matrix element. 

In Fig. [5] the differences between the CCM and the 
Miller-Spencer SRC are manifested by plotting the ratio 
of correlated H src+ f ns (r 12) and uncorrelated Hb are (r 12) 
neutrino potentials. The averaged energy of intermediate 
nuclear states E is 8 MeV. For sake of simplicity the effect 
of higher order terms of nucleon currents on the neutrino 
potential is neglected. From Fig. [2] we see that there is 
a significant difference between the CCM and the Miller- 
Spencer treatment of the SRC. The maxima of the CCM 
and the Spencer-Miller curves occur at 1 fm and 1.5 fm, 
respectively. One finds also that the reduction at short 
distances is much weaker for CD-Bonn than for Argonne 
interactions. 

For purpose of numerical calculation of the 0^,3/3-decay 
NMEs we present the CCM short-range correlation func- 
tions in an analytic form of Jastrow-like function as 



fAMm) = 1 - ce- ar -(l-br 2 ). 



(22) 
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FIG. 3: The r i2 dependence of M°" in 76 Ge evaluated in the 
model space that contains twelve subshells. The four curves 
show the effects of different treatment of short-range corre- 
lations (SRC). The finite nucleon size (FNS) is taken into 
account. 



C. Finite nucleon size and two-body short range 
correlations 



The FNS effects are introduced in the calculation of 
the Ov (3 /3-decay NMEs by the dipole form factors in mo- 
mentum space. The form factor simulates the fact that 
nucleon is not a point particle, and therefore as q in- 
creases, the probability that nucleon will stay intact (and 
not produce pions etc) decreases. The physics of FNS 
and SRC is different, but both reduce the magnitude of 
the operator when q 2 increases or equivalently r\ 2 de- 
creases. It was found [2(| that Miller-Spencer and the 
UCOM short-range correlations essentially eliminate the 
effect of the FNS on the matrix elements. The same is 
expected to be valid also for the CCM CD-Bonn and Ar- 
gonne short-range correlations. From Fig. [5] we see that 
the ratio of correlated and uncorrelated neutrino poten- 
tials is changed only weakly if in addition to two-nucleon 
SRC the effect of the FNS is taken into account. 



The set of parameters for Argonne and CD-Bonn NN 
interactions is given by 

Ja{txi): a = 1.59 fm- 2 , b = 1.45 /m~ 2 , c = 0.92, 
/s(ri 2 ) : a = 1.52 fm- 2 , b = 1.88 fm-' 2 , c = 0.46. 

(23) 

The calculated NMEs with these short-range correlation 
functions agree within a few percentages with those ob- 
tained without this approximation. We note that the 
dependence of the SRC on the value of oscillator length 
b is rather weak. 

In Fig. [3]the r\ 2 dependence of M 0l/ is shown for CCM 
Argonne, CCM CD-Bonn and phenomenological Jastrow 
SRC for the O^/3/3-decay of 76 Ge. The quantity C{r 12 ) is 
defined by 

/>oo 

M 0v = / C(r 12 )dr 12 . (24) 
Jo 

We note that the range of r\ 2 is practically restricted 
from above by r\ 2 < 2R nuc i. From the Fig. [3] we see that 
a modification of the neutrino potential due to the differ- 
ent types of SRC is transmitted to the behavior of C(ri 2 ) 
for ri2 < 2fm. Both, the CCM short-range correlation 
functions (see Fig. [2]) and C(r 12 ) with SRC switched 
off (but with FNS effect) have maxima for r\ 2 ~ 1 fm 
unlike the phenomenological Jastrow function with max- 
imum shifted to r\ 2 = 1.5 fm. This explains a signifi- 
cant increase of C(r\ 2 ) with CCM SRC and suppression 
of C(ri 2 ) with Jastrow SRC in this region. This phe- 
nomenon clarifies also why the values of M 0u obtained 
with CCM SRC are comparable to those calculated when 
only the FNS effect is considered (see Table U). The in- 
crease of C(ri 2 ) for ri2 ~ 1 fm compensates its reduction 
in the range r\ 2 < 0.7 fm. 



It is worth mentioning that the behavior of the UCOM 
correlated neutrino potential differs strongly from those 
calculated with the CCM and Jastrow SRC. This is man- 
ifested in Fig. 2J The studied ratio of UCOM correlated 
and uncorrelated neutrino potentials never exceeds unity 
unlike in the case of CCM correlations (see Fig. fj]). Ac- 
tually, the UCOM SRC imitate the effect of the FNS with 
form-factor cut-off of about 850 MeV. The two-nucleon 
wave function can be treated as two point-like objects for 
nucleon separations greater than about 1.5 fm. 

The effect of the SRC on the 0^^/3-decay NMEs has 
been referred mostly to case when the FNS is taken into 
account. It was found that Miller-Spencer SRC reduces 
the 0*W35-decay NMEs by 20-30% and UCOM SRC by 
~5% [201 ■ F° r better understanding of this effect we cal- 
culate the 0^/3-decay of 76 Ge, 1Q0 Mo and 130 Te with 
and without consideration of the FNS. The 12 levels 
( 76 Ge) and 13 levels ( 100 Mo and 130 Te) single particle 
model spaces are used in calculation. The results are 
displayed in Table |U The bare NME was obtained in 
the limit of cut-off masses My, a go to infinity and with 
the two-nucleon SRC switched off. The FNS values of 
M 0v are determined by nucleon form factors with phe- 
nomenological values of My and Ma- We see that the 
FNS reduces M 0v by 20 %. The O^/3/3-decay NME is 
suppressed by about 20-30% and 40-50% in the cases 
of the two-nucleon CCM (Argonne potential) and the 
phenomenological Miller-Spencer SRC, respectively. It is 
also shown that once SRC effects are included the con- 
sideration of the nucleon form-factors almost does not in- 
fluence the value of M ". It is because the FNS and the 
SRC effects act coherently on the Ov ft (3-decay NMEs and 
both diminish them. However, the effect of the SRC is at 
least partially stronger (CCM SRC) or stronger (Miller- 
Spencer SRC) in comparison with the FNS effect. 
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1.2 ^ 




r 12 [fm] 

FIG. 4: A ratio of neutrino potentials with inclusion of the 
finite nucleon size effect (FNS) and bare neutrino potential 
for cut-off masses My and My . This is compared with ratio 
of neutrino potential with UCOM short range correlations 
(SRC) and bare neutrino potential. E = 8 MeV is used in 
calculation. 



IV. NUMERICAL RESULTS 

The nuclear matrix elements of the O^/3/3-decay of the 
experimentally interesting nuclei A — 76, 82, 96, 100, 
116, 128, 130 and 136 are systematically evaluated using 
the QRPA and RQRPA. In the present calculations, we 
improve on the Miller-Spencer Jastrow and UCOM meth- 
ods by engaging the SRC calculated within the exp(S) 
approach with the CD-Bonn and Argonne V18 NN inter- 
actions. This allows a consistent study of the Oz^/3/3-decay 
NMEs as for the first time the same realistic nucleon- 
nucleon force is used for the description of the pairing 
interactions, RPA ground state correlations and the two- 
nuclcon SRC. 

The nuclear structure calculations are p erformed as de- 
scribed in our previous publications [la, [l^, [2(| • Three 
different single-particle model spaces are used: small (2- 
3 oscillator shells), intermediate (3-4 oscillator shells) 
and large (5 oscillator shells) model spaces (see Ref. 
[ll|). The single-particle energies are obtained by using a 
Coulomb-corrected Woods-Saxon potential [4l[. The in- 
teractions employed are the Brueckner G-matrices which 
are a solution of the Bethe-Goldstone equation with the 
CD-Bonn and Argonne V18 one-boson exchange poten- 
tials. The pairing two-body interaction is fitted in the 
standard way and the pairing parameters of the BCS 
are adjusted to reproduce the phenomenological pairing 
gaps, extracted from the atomic mass tables. We renor- 
malize the particle-particle and particle-hole channels of 
the G-matrix interaction of the nuclear Hamiltonian by 
introducing the parameters g pp and g p h, respectively. We 
use g p h = 1 throughout what allows to reproduce well the 
available data on the position of the giant Gamow- Teller 
resonance. The particle-particle strength parameter g pp 
of the (R)QRPA is fixed by the data on the two-neutrino 



double beta decays. 

The NME calculated within the above procedure, 
which includes three different model spaces, is denoted as 
the averaged 0z//3/?-decay NME (M' 0l/ ). The results are 
presented separately for the CD-Bonn and Argonne inter- 
actions and for two different values of the axial coupling 
constant gA in Table [Til We confirm again that with con- 
sidered procedure the O^/3/3-decay values become essen- 
tially independent on the size of the single-particle basis 
and rather stable with respect to the possible quenching 
of the g A . The NMEs obtained with the CD-Bonn NN 
interaction are slightly larger as those calculated with the 
Argonne interaction. This is explained by the fact that 
the CCM Argonne correlation function cuts out more the 
small ri2 part from the relative wave function of the two- 
nucleons as the CCM CD-Bonn correlation function. The 
differences in NMEs due to a different treatment of the 
SRC do not exceed differences between the QRPA and 
the RQRPA results. 

In Table IIIII we show the calculated ranges of the nu- 
clear matrix element M 0v evaluated within the QRPA 
and RQRPA, with standard (gA — 1.254) and quenched 
(gA = 1.0) axial- vector couplings and with the CCM CD- 
Bonn and Argonne SRC functions. These ranges quantify 
the uncertainty in the calculated O^/3/3-decay NMEs. In 
respect to the central value their accuracy is of the order 
of 25%. A significant amount of the uncertainty is due to 
a quenching of the axial-vector coupling constant g A in 
nuclear medium [42| . For a comparison we present also 
the NMEs calculated with the phenomenological Jastrow 
SRC function in Table IIIII The notable differences be- 
tween the results calculated with Jastrow and CCM SRC 
are about of 20-30%. Of course, the results obtained with 
the CCM SRC are preferable. Unfortunately, they can 
not be directly compared to those of the complementary 
Large-Scale Shell Model (LSSM) as they have been eval- 
uated only with the Jastrow and UCOM SRC and for 
gA = 1-25 [HI, Hlj]. It is reasonable to assume that the 
LSSM values (see Table UTTJ would be increased also by 
the order of 20-30%, if the CCM SRC would be consid- 
ered. Recall, that within both approaches qualitatively 
the same ri2 dependence of M 0v was found (20.[2l|. Fig. 
[5] shows our calculated ranges for M' 0l/ , which are com- 
pared with ranges calculated with Miller-Spencer Jastrow 
function and the latest LSSM results. Given the interest 
in the subject, we show the range of predicted half-lives 
corresponding to our full range of M 0v in Table IIIII 

Recently, the occupation numbers of neutron and pro- 
ton valence orbits in the 76 Ge and 76 Se nuclei were mea- 
sured by neutron and proton adding and removing trans- 
fer reactions [HI, • In following theoretical study (45[ 
these results were used as a guideline for modification 
of the effective mean field energies that results in bet- 
ter description of these quantities. The calculation of 
the O^/3/3-decay NME for 76 Ge performed with adjusted 
mean field [45| and in combination with the selfconsistent 
RQRPA (SRQRPA) method [H, which conserves the 
mean particle number in correlated ground state, lead to 
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FIG. 5: The full ranges of M' " with the CCM and Miller- 
Spencer Jastrow treatments of the short range correlations. 
For comparison the results of a recent Large Scale Shell Model 
evaluation of M " that used the Jastrow-type treatment of 
short range correlations are also shown. 



a reduction of M by 25% when compared to previous 
QRPA value. The phenomelogical Jastrow and UCOM 
SRC were considered. We found that a reduction of 20% 
in the case of CCM CD-Bonn SRC. We have 



< M 0l/ > 



= 4.24(0.44), 3.49(0.23), Jastrow 
= 5.19(0.54), 4.60(0.23), UCOM 
= 6.32(0.32), 5.15(0.44), CCM. 



the 0i//3/3-decay. These correlations, which have origin 
in the short-range repulsion of the realistic NN inter- 
action, are neglected in the mean-field, the LSSM and 
the QRPA descriptions. Till now, Miller-Spencer Jas- 
trow and the UCOM SRC have been introduced into the 
involved two-body transition matrix elements, changing 
two neutrons into two protons, to achieve healing of the 
correlated wave functions. The effect of these SRC was 
considered as an uncertainty [20j. 

In this article the CCM short-range correlations has 
been considered. They were obtained as a solution of 
the coupled cluster method with realistic CD-Bonn and 
Argonne V18 interactions. An analysis of the squared 
correlation functions, represented by a ratio of correlated 
and uncorrelated neutrino potentials, has showed a prin- 
cipal differences among the Miller-Spencer, UCOM and 
CCM SRC. In addition, the importance of the effect of 
the FNS was studied. It was found that both CCM SRC 
and the FNS effect reduces the Oi//3/3-decay NMEs by 
a comparable amount for a considered choice of form- 
factor masses. The suppression due to Miller-Spencer 
SRC is about twice larger when compared to results with- 
out SRC and the FNS effect. 

Finally, we have improved the presently available cal- 
culations by performing a consistent calculation of the 
Oi//3/3-decay NMEs in which pairing, ground-state corre- 
lations and the short-range correlations originate from 
the same realistic NN interaction, namely from the CD- 
Bonn and Argonne potentials. 



(25) 



Here, the first and second values in each line correspond 
to the QRPA with Woods-Saxon mean field and the 
SRQRPA with adjusted Woods-Saxon mean field [Ifl 
way of calculations, respectively. In order to better un- 
derstand the role of proton and neutron occupation num- 
bers in the Oi//3/3-decay calculation further experimental 
and theoretical studies are needed what goes beyond the 
scope of this paper. 



V. CONCLUSIONS 

We have addressed the issue of a consistent treat- 
ment of the short-range correlations in the context of 
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TABLE I: Nuclear matrix elements for the O^/3/3-decays of 76 Ge, 100 Mo and 136 Te within the QRPA. The results are presented: 
i) (bare) no correlations and no nucleon form factors; ii) (FNS) no correlations but with nucleon form factors; iii) (SRC) CCM 
Argonne and Miller-Spencer short-range correlations but without nucleon form factors; iv) (FNS+SRC) correlations and nucleon 
form factors. 



Nucleus 


bare 


FNS 


SRC 


FNS + SRC 








CCM Miller-Spencer 


CCM Miller-Spencer 


ioo Mo _> ioo Ru 

130 Te ^ 130 Xe 


7.39 
6.15 
5.62 


6.14 
4.75 
4.49 


5.86 4.46 
4.40 2.87 
4.22 2.97 


5.91 4.54 
4.46 2.96 
4.27 3.04 



TABLE II: Averaged 0v{3(3 nuclear matrix elements {M' v ) and their variance a (in parenthesses) calculated within the RQRPA 
and the QRPA. The pairing and residual interactions of the nuclear Hamiltonian and the two-nucleon short-range correlations 
(SRC) are derived from the same realistic nucleon-nucleon interaction (CD-Bonn and Argonne potentials) by exploiting the 
Brueckner-Hartree-Fock and coupled cluster (CCM) methods. Three sets of single particle level schemes are used, ranging in 
size from 9 to 23 orbits. The strength of the particle-particle interaction is adjusted so the experimental value of the 2v(5[3- 
decay nuclear matrix element M^f is correctly reproduced. Both free nucleon (gA = 1.254) and quenched (gA = 1-0) values of 
axial- vector coupling constant are considered. 



Nuclear 
transition 


9A 


[MeV" 1 ] 


(M' u,y ) 


{M' { " 





CCM CD-Bonn SRC 


CCM Argonne SRC 


RQRPA 


QRPA 


RQRPA 


QRPA 


ib Ge -> ' ,b Se 


1.25 


0.15 


5.44(0.23) 


6.32(0.32) 


4.97(0.19) 


5.81(0.27) 




1.00 


0.23 


4.62(0.22) 


5.16(0.25) 


4.21(0.14) 


4.77(0.20) 


M S e - 82 K r 


1.25 


0.10 


4.86(0.20) 


5.65(0.27) 


4.44(0.19) 


5.19(0.24) 




1.00 


0.16 


3.93(0.15) 


4.48(0.20) 


3.67(0.14) 


4.19(0.18) 


^Zr -> 96 Mo 


1.25 


0.11 


2.01(0.20) 


2.09(0.03) 


1.84(0.16) 


1.90(0.09) 




1.00 


0.17 


1.72(0.15) 


1.93(0.11) 


1.55(0.12) 


1.74(0.11) 


ioo Mo ^ ioo Ru 


1.25 


0.22 


4.28(0.28) 


5.25(0.31) 


3.85(0.31) 


4.75(0.33) 




1.00 


0.34 


3.44(0.19) 


4.07(0.22) 


3.14(0.23) 


3.69(0.25) 


^Cd~, 116 Sn 


1.25 


0.12 


3.41(0.24) 


3.99(0.15) 


3.06(0.22) 


3.54(0.27) 


128 Te 128 Xe 


1.00 


0.19 


2.68(0.19) 


3.03(0.19) 


2.47(0.17) 


2.74(0.21) 


1.25 


0.034 


4.82(0.15) 


5.49(0.16) 


4.32(0.16) 


4.93(0.16) 


130 Te 130 Xe 


1.00 


0.053 


3.67(0.11) 


4.16(0.12) 


3.32(0.11) 


3.77(0.12) 


1.25 


0.036 


4.40(0.13) 


4.92(0.12) 


3.91(0.14) 


4.37(0.14) 




1.00 


0.056 


3.38(0.08) 


3.77(0.07) 


3.02(0.10) 


3.38(0.10) 


136 Xe ^ 136 Bfl 


1.25 


0.030 


2.89(0.17) 


3.11(0.13) 


2.59(0.16) 


2.78(0.13) 




1.00 


0.045 


2.26(0.11) 


2.42(0.08) 


2.03(0.10) 


2.17(0.09) 




1.25 





2.53(0.17) 


2.73(0.13) 


2.25(0.16) 


2.43(0.13) 




1.00 





1.87(0.11) 


2.01(0.08) 


1.67(0.10) 


1.80(0.09) 



11 



TABLE III: The calculated ranges of the nuclear matrix element M evaluated within the QRPA and RQRPA, with standard 
(gA = 1.254) and quenched (qa = 1.0) axial-vector couplings and with coupled cluster method (CCM) CD-Bonn and Argonne 
short-range correlation (SRC) functions. Column 4 contains the ranges of M °" with the phenomenological Miller-Spencer 
Jastrow treatment of short range correlations, while column 6 shows the CCM SRC-based results. For comparison the recent 
results of a large scale shell model (LSSM) evaluation of M' 0l/ [l|| that used the Miller-Spencer Jastrow SRC and gA = 1-25 
are given in column 2. However, they have to be scaled by factor (1.1 /m/1.2 fm) as different value of ro (Rnuci = tqA 1 ^ 3 ) was 
considered. Columns 3, 5 and 7 give the 0^/3/3-decay half-lives or half-life ranges corresponding to values of the matrix-elements 
in columns 2, 4 and 6 for < mpp >= 50 meV. 



Nucleus LSSM (Jastrow SRC) (R)QRPA (Jastrow SRC) (R)QRPA (CCM SRC) 





M 0v T°J 2 


{(mpp) 


= 50 meV) 


M'°" 


T?J 2 {(mpp) = 


: 50 meV) 




T?/ 2 ({mpp) = 


: 50 meV) 


' lti Ge 


2.22 


3.18 x 


10 2V 


(3.33,4.68) 


(6.01, 11.9) 


x 10^ b 


(4.07,6.64) 


(2.99, 7.95) 


x 10™ 


S2 Se 


2.11 


7.93 x 


10 26 


(2.82,4.17) 


(1.71,3.73) 


x 10 26 


(3.53,5.92) 


(0.85,2.38) 


x 10 26 


96 Zr 








(1.01, 1.34) 


(7.90, 13.9) 


x 10 26 


(1.43,2.12) 


(3.16,6.94) 


x 10 26 


100 Mo 








(2.22,3.53) 


(1.46,3.70) 


x 10 26 


(2.91,5.56) 


(0.59,2.15) 


x 10 26 


116 Cd 








(1.83,2.93) 


(1.95,5.01) 


x 10 26 


(2.30,4.14) 


(0.98,3.17) 


x 10 26 


i28 Te 


2.26 


1.10 x 


10 28 


(2.46, 3.77) 


(3.33, 7.81) 


x 10 27 


(3.21,5.65) 


(1.48,4.59) 


x 10 27 


i30 Te 


2.04 


5.39 x 


10 26 


(2.27,3.38) 


(1.65,3.66) 


x 10 26 


(2.92,5.04) 


(7.42,2.21) 


x 10 26 


136 Xe 


1.70 


6.79 x 


10 26 


(1.17,2.22) 


(3.59, 12.9) 


x 10 26 


(1.57,3.24) 


(1.68, 7.17) 


x 10 26 



